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GENERAL NOMENCLATURE 


I 
- dynamic pressure impulse 


D * displacement 


R 
ground range 


x a scaled ground range 


HOB m height of burst 


Sd , Sit Sp, St a scaling factors 


Po 
1/3 
1 
sd -(T4 


T 
.(*Q 
+ 273 1/2 
14,7 2/3 


14.7 
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00 
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14.7" 


To(C) * 273 1/2 
Po 
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11/3 
st % 2-,33 
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P0 a ambient pressure (lbs/sq in) 
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ambient temperature (degrees C) 


V a weapon yield (KT) 


aO, a 1 , a2, b, t, a, A, B a constants in least squares fits 


x,y - independent and dependent variables in least squares fits 


a standard deviation 


a 
variance 


ai 


N * number of data points 


r a correlation coefficient 


R • multiple linear correlation coefficient 
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GENERAL NOMsENCLATURE (Continued) 


1.,R coefficient of determination 


S *sun 
of squares of deviations of data points from fitted curve 


S' *sum of squares of fractional (relative) deviations o~f data points from 


fitted curve 


nubo 
1/2 w standard deviation of 
(N4 
- nubrof regression coefficieonts 
fractional error 


St1/2 
I * -i-) 
x 100 
'~root 
mean square percent error 


SO - side on 


FO a face an 


NO a rear on 


See Appendix, "The Least Squares Procedures", for definitions/defining equa- 


tions of the following specific examples of the above general quantities: 
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SECTION I 


(This Section is Unclassified) 


INTRODUCTION 


The objective was to construct improved height of burst contours for 


dynamic pressure impulse. 
In particular, we wished to extend the set of con- 


tours previously generated (Reference 1) to higher values of dynamic pressure 
impulse and to obtain data points along the contours at additional scaled 
burst heights, i.e., to better determine the contour shapes. To do this we 
analyzed all of the data available. This includes dynamic pressure-time wave- 
forms for shots not used in Reference I (because no trucks were exposed on these 
shots): TEAPOT Hornet-S, Post-ll, and Zucchini-14; PLIMBBOB Franklin-2, Wilson-4, 
Hood-6, Kepler-9, and Owens-tO. It also includes data for shots in which tanks 
and howitzers were exposed but no dynamic pressure measurements were made: 
GREENHOUSE Easy-2; TUMBLER-SNAPPER Fox-6 and How-?; UPSHOT-INOTHOLE Annie-l*, 
lp 


Nancy-2, Badger-S, and Simon-?. 


-• .- 


I 


There was one dynamic pressure measurement on this shot. 
, 
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SECTION 2 


(This Section is Unclassified) 


ANALYSIS OF BLAST MEASUREMENTS 


2.1 
Measurements For Evtv.tp Not Previously Considered. 


We consider first the shots not previously considered for which some 


dynamic pressure measurements were available. Some of the waveforms were not 
very good owing to noise (mixed inextricably with real physical irregularities), 
baseline movement, or gage record cutoff. However, many of the waveforms were 
found useable, being a quality similar to, or not significantly inferior to, 
those used in Rr.ference 1. General information on all of the nuclear tests of 
concern in this report is given in Table 1. Dynamic pressure impulse, scaled 
and unscaled, is given in Table 2 for each Operation/Event, scaled height of 
burst, and ground range for which we have obtained results. Also given in 
Table 2 is the manner in which we obtained the values of dynamic pressure 
impulse. 


We discuss the data for the shots in the order listed in Section 1. 


TEAPOT, Hornet-5. As indicated in Table 1 we used BRL gage data.* 


There were two gage results at each of the three ground ranges. 
In two instan- 


ces the two results are in good agreement. At the closest-in range, 256 meters, 
the two gages differ by about a factor of two but one of the two appears to be 
correct, the other badly in error. We used the average of the two results for 
the 329 and 460 metre ranges, and the apparently correct results for the 256 
metre ground range. 


TEAPOT, Post-ll. For this shot dynamic pressure waveforms are avail- 


able in Reference 2. We used a planimeter to obtain the area under each curve 
and then used the data reduction procedure described in detail in Reference 1. 
A detailed discussion of accuracy of this procedure is also given in Reference 
1. For this shot waveforms are given at four ground ranges in Reference 2. 
The waveforms are not good but are of sufficient quality to be useful. 
As 


indicated in Table 2 we used the average of BRL's result and our result for 
the two ground ranges at which BRL provided results. 
Our resulto and BRL's 


results agreed to within about 3% which is quite good; (the two results are 
based on the same gage data but the data reduction procedures are indepeident; 
we have analyzed errors inherent in data reduction in Reference 1, and, in 
many instances, errors can be quite large). 


For the two ground ranges at which BRL did not 


reduce the data, we used the result of our own data reduction process. 


TEAPOT, Zucchini-14. For this shot we averaged our result with 


the result of BRL's data reduction process for the three listed ground ranges. 


* BRL gage data used in this report are published In DNA-TR-85-161. 


Furnished by Mr. J. Keefer, BRL. 


!K 


8 


UNCLASSIFIED 


) 
%7= 


UNCLASSIFIED 


In each case the difference was between 3 and S%. 
Results are based on two 


gages at the 610 and 700 metro ground ranges and upon one gage at the 794 
metre range. 


PLUMBBOB, Franklin-2. For this shot we used the BRL result for the 


single ground range for which there appeared to be a good measurement. 


PLU BBOB, Wilson-4, For this shot BRL also provided values of 


dynamic pressure impulse as listed in Table 2. (We integrated the waveform 
curves as a check, however.) 


PLUMBBOB, Hood-6. Of the 6 ranges at which there were gage results 


the data appeared to be valid at 4 locations. There were two gages at each of 
these locations. We used the BRL results as listed in Table 2 (again after an 
integration check). 


PL14BBOB, 
Kepler-9. In this case waveforas obtained at three ground 


ranges were quite poor. In one case there was a large disparity between two 
BRL gage results (almost a factor of 4); in a second case there was almost a 
factor of 2 difference between our gage-reduced result and BRL's result. 
At 


the third ground range, 762 metres, the waveform was somewhat better. 
The 


percent difference between our result and BRL's result was rather large, about 
20%; in Table 2 we list the average of the two results. We discarded the 
results at the other two ground ranges. 


PLUMOBOB, 
Owens-lO. 
For this shot we again used the average of our 


gage-reduced result and the BRL result at two of the three ground ranges, 305 
and 518 metres. At the 305 metre range the results are in reasonably good 
agreement (within 9%); at the 518 metre range the two results differ by a 
factor of 1.4 but the waveforms are quite poor. At the 457 metre range our 
value of dynamic pressure impulse was exactly equal to BRL's value of scaled 
dynamic pressure impulse; the BRL value appeared anomalous on a data plot, 
probably as a result of failure to scale the value obtained. In this case we 
used our gage-reduced result. 


2.2 Analysis of Tank Data. 


We next consider the set of shots listed in Section 1 on which tanks 


were exposed but there were no dynamic pressure measurements. Our rationale 
ell 


here is based on a finding we have discussed in detail in Reference 1: 
the 
1P 


displacement which a vehicle exposed to a blast wave sufers can be used as 
a measure of the dynamic pressure impulse it receives. 
,hat is, the vehicle 


can serve as a gage for dynamic pressure impulse if we are able to "calibrate" 
this "gage". The calibration is the curve of dynamic pressure impulse versus 
displacement. 
(We are ignoring low yield devices where diffraction effects 


play a role.) 
We have analyzed the procedure and results for 1/4 ton and 2 


ton trucks in considerable detail in Reference 1. 


We would not be able to use the procedure for tanks if we had no dynamic 


pressure impulse data for the ground ranges at which tanks were exposed. That 
is, we would have no way of obtaining the necessary calibration curve. For- 
tunately, we have the necessary dynamic pressure impulse data. It is supplied 
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by information from TEAPOT, Met-12 and Apple 11-13 and from UPSHOT-KNOTHOLE, 
Annie-i and Grable-lO. The data for the pair of UPSHOT-MOTHOLE events are 
gage data supplied by BRL; for Grable-lO we also have results based on both 
gage and (truck) displacement from Reference 1, but at larger values of ground 
range than those at which tanks were exposed. The BRL data and the Reference 1 
data are very compatible, i.e., very well represented by a single curve of 
scaled dynamic pressure impulse versus scaled ground range (as shown in Figure 
7, to be discussed later). The results are listed in Table 2. 


For Met-12 we have the results from Reference 1, listed in Table 2, based 


on both gage and (truck) displacement data. (There is a considerably larger 
body of truck displacement data than of tank displacement data so that the 
dynamic pressure impulse-displacement curve, i.e., the calibration curve, is 
better determined for trucks.) 


For Apple 11-13 the results listed in Table 2 are taken only partly from 


Reference 1. In this instance, values of dynamic pressure impulse at the 518 
and 625 metre ranges listed in Reference 1 (based only upon gage data) have 
been averaged with BRL gage results, i.e., the two results given equal weight. 
At the 625 metre range this.makes a negligible difference while at the S18 
metre range the value of dynamic pressure impulse listed in Reference 1 differs 
fram that of Table 2 by about 16%. Finally at the 808 metre range the result 
listed in Table 2 differs from that listed in Reference 1; the latter is in 
error as a check of our (previous) data reduction reveals. 


2.2.1 Dynamic Pressure Impulse Versus Displacement For Tanks. 


We plotted dynamic pressure impulse versus displacement for the 


tanks exposed on events Met-12, Apple 11-13, Annie-1 and Grable-10. This 
provides a total of 14 data points* listed in Table 3 along with the tank 
displacement data for other shots. 


There are several variables which cause deviations of the data 


relative to a smooth curve through the points: 
(1) vehicle orientation with 


respect to the bomb -- 4 data points correspond to side-on orientation, 7 to 
face-on orientation, 1 to rear-on orientation, 1 to face-on 45 orientation, 
and 1 to face-on 3/4 left orientation; (2) there are three different tanks -- 
M4A3, M24, and M48; (3) there are two surface conditions -- rough sand and 
fine sand. There obviously are not sufficient data to sort out the effects 
of the variables. However, when we plotted the data we found that: 
(1) 


(initial) orientation appears to make little difference, any effect being 
submerged in effects of the other variables; (we intuitively expect initial 
T.,;1 


orientation to make less difference for tanks than for trucks); (2) any 
systematic deviation which could be attributed to type of vehicle could as 
well be attributed to inaccuracy in dynamic pressure impulse -- for example, 
. 
L 


the Met-12 points (points B, 6, and 17 in Table 3) are a little high and are 


- I 


-We havenot used Point 8 of Table 3, a Smoky-IS point. Placing this point 
on Figure 1 shows that there is clearly something wrong with it -- and what 
is wrong with it involves the displacement, not the dynamic pressure impulse, 
', 


even though the latter involves extrapolation on Figure 6. 
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M48 points while the Apple 11-13 points are a little low and pertain to both 
M48's and 24's. 
(In fact, if we use an average dynamic pressure impulse 


versus scaled ground range curve for the two shots, Met-12 and Apple 11-13, 
these deviations disappear.) 


Plots of the data and the fits obtained are shown in Figure i.* 


We used these data to infer dynamic pressure impulse from the 


measured displacements on events GREENHOUSE Easy-2; UPSHOT-KNOTHOLE Nancy-2, 
Badger-5, and Simon-7; TUMBLER-SNAPPER Fox-6 and How-7. 


We have noted in Reference I that, despite rather wide fluctuations 


in displacement owing to the several variable factors mentioned, the central 
curve following the trend of the data points provides a rather good calibration 
which enables us to use a measured displacement to determine the dynamic pres- 
sure impulse to which the vehicle was subjected. 
Since displacements vary owing 


to uncontrolled factors, the inferred dynamic pressure impulses will exhibit 
f 


considerable dispersion. 
(However, in many cases described in Reference 1, the 


displacement-inferred values are about as reliable as the gage-inferred values, 
especially when there are two or more values which can be averaged.) 


The calibration curve for tanks is less well determined than that 


for trucks (Reference 1), there being fewer data points. Also the displacements 
are smaller resulting in greater errors especially from small yield. Nonetheless, 
when the displacement-inferred dynamic pressure impulse values are plotted versus 
scaled ground range (along with a few deta peints for which there also are gage 
data), we see that the curve obtained is reasonably well determined. 
See Figure 


3. 


We used an eye-drawn curve rather than a proportional fit in Figure 1 


in inferring dynamic pressure impulse from measured displacements for the follow- 
ing reasons: 
J 


(1) point 20, very small displacement, is much less important 


than the other data points in application of this fit, i.e., inference of 
dynamic pressure impulse from measured displacement; 


(2) point 21 has virtually the same displacement as the cen- 


troid so that it can be ignored in drawing a straight line through the data 
points; [the line must pass through the centroid: 


r-IT= 2.448, r-T- = 3.768; 
" 


q 


on the abscissa and ordinate scales from Figure 1 we see that both fits do 
pass through the point D(nD) - 11.56, Iq(Zn lq) = 43.29; the lines do not 


q 


pass through the point 
= 21.85, I-= 47.95]; 


q 


* Because of the mentioned variables whose effects we cannot disentangle owing 


to the small amount of data, we have numbered the points on Figure 1. Cor- 
responding numbers are listed in Table 3. By comparison one can verify 
that the data do not exhibit any marked effects owing to differences among 
the above variables. 
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(3) the Grable-lO points 21 and 22 are badly inconsistent with 


one another and points 9 and 10 are somewhat inconsistent. Thus our eye-diawn 
curve was a little higher than Fit 1 of Figure 1 at the small displacement end 
and, like Fit 1, passed through the centroid of the data points. The differ- 
ence, however, is nowhere greater than a few percent. We also note from 
Figure 1 that a curve through the data points would do no better than a 
straight line and is therefore not warranted. It is quite possible, however, 
that if the data covered a wider range a straight line would no longer be 
adequate. 


2.3 
Discussion Of and Minor Improvements Upon Previously Used Blast Data. 


Most of the dynamic pressure impulse data listed in Table 2 which we 


have not yet discussed is taken from Reference 1. In a few instances a value 
taken from Reference 1 (gage only - no vehicle displacement data available) 
was averaged with a BRL gage-reduced result as shown in the table. 
In one 


instance, Turk-4, ground range 59S metres, we reduced the gage data even 
though the gage had cut off near the end. We estimated the shape and time of 
pulse completion from other Turk-4 data and believe that the inaccuracy so 
introduced is small - 
not more than a few percent. Reference to Figure 3 


shows why this is desirable: Turk-4 has the lowest scaled burst height for 
which we have gage data (except for a single Annie-1 point and except for 
surface burst data); its scaled burst height is only moderately greater than 
those of other events shown on Figure 3; the Turk-4 data are consistent with 
the data for the other events shown on Figure 3 and extend to a much greater 
value of scaled ground range. 
Thus the curve fit of Figure 3 is not diminished 


in reliability by inclusion of Turk-4 data which, however, enable its use over 
a much greater range in constructing height of burst charts, the object of this 
report. 


For events Yuma-4 and Wasp Prime-9 we used our values from Reference 1, 


therein described as "first iteration" results; i.e., when available, truck 
displacement-inferred dynamic pressure impulse values were averaged with gage 
values in obtaining our best estimates of dynamic pressure impulse. 
(The 


manner of averaging, justification, and discussion of accuracy and reliability 
are considered in detail in Reference 1.) 
We originally considered Yuma-4 to 


be an ideal event, but we now believe the values achieved by averaging displace- 
ment-inferred dynamic pressure impulse values with the gage values represent an 
improvement over the gage values alone; the results achieved are certainly more 
compatible with the other data in the same scaled height of burst region. See 
Figure 9 (and compare with gage only data listed in Reference 1, page 11S). 


For event Encore-9, a near-ideal event, we use the gage results -- even 


though in Reference 1 we sought to improve upon these results by use of 
displacement-inferred dynamic pressure impulse. The reasoning used in Reference 
1 was that since the data scatter on a displacement versus dynamic pressure 
impulse plot was just as great for ideal/near-ideal shots as for non-ideal and 
since the displacement-inferred dynamic pressure impulse for non-ideal shots 
is, on the average, as accurate and reliable as the gage values -- then 
averaging of displacement-inferred with gage-inferred dynamic pressure impulse 
should improve the values for near-ideal events just as it does for non-ideal. 
The points at issue are discussed in detail in Reference 1. One point, however, 
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is that there is a much greater variety of ideal/near-ideal shots than of 
non-ideal, especially regarding the range of weapon yields in the data base.' 
When good gage data are available, these data should produce better results 
than can be obtained by averaging in displacement-inferred values -- owing to 
the effects of the several uncontrolled variablez previously described. That 
is, while, on the average, displacement-inferred values are reliable, consider- 
able dispersion is to be expected. Compare, for example, Figure 3; here with 
the exception of the Turk-4 data points and a single Annie-1 point, all of the 
plotted points are based entirely upon displacement-inferred values of dynamic 
pressure impulse. While the central curve following the general trend of the 
data is reasonably well determined, the data point deviations are rather large 
N. 


-- 
considerably larger than the deviations in Figures 4-10. 


Near-ideal shots also have (small) real differences relative to one 


another. This is shown, for example, by"Figure 3.13 of Reference 1 which is 
a plot of scaled dynamic pressure impulse versus scaled ground range. Data 
for the near-ideal (surface burst) shots do not completely coalesce under the 
scaling -- there are clearly small but real differences owing to variable 
factors not controlled in the experiments. In Figures 3-10 the data for the 
events plotted on each figure seem to coalesce quite well with respect to a 
single curve. In Figure 7 some deviation can be seen. In general, deviations 
tend to appear when data for several events are shown on a single plot, when 
there are several data points for each event, and where the data for several 
events overlap, i.e., cover the same domain of the abscissa, scaled ground 
range. The data in Figures 3-10 per:ain to a much smaller range of weapon 
yields than do the surface burst dati of Reference 1 (Figures 3.13 and 3.15). 
% 


In general, when the data for sever; & shots coalesce we can attach a high 
degree of reliability to the data. 


.J4i 


In constructing height of burst charts in this report we used the results 


described and listed in Table 3. For the surface burst data we used Figure 
3.13 of Reference 1. 


Finally we analyzed field data for many more shots, but found the data 


unusable. In some instances the waveforms were very bad; in other instances 
the gages were placed at elevations other than 3 feet above the ground. The 
effect of gage elevation is non-negligible.- (Different types of gages also 
exhibit somewhat different responses which must be accounted for in order to 
achieve consistent results.) 


NOTE: The BRL data used in Table 3 have not been published. 


2.4 Dynamic Pressure Impulse Versus Displacement for Self-Propelled Howitzers. 


Figure 2 is a plot of displacement - dynamic pressure impulse data for 


For a single (extrapolated) point on the charts for which the scaled dynamic 
pressure impulse is IS kPa-sec and the sealed height of burst "is zero, we 
also used Figure 3.15 of Reference 1 and took account of the difference 
between Figures 3.13 and 3.15 in this region. Figure 3.15 extends to slight- 
ly higher values of dynamic pressure impulse than does Figure 3.13 so that 
less extrapolation is required. 
-: 


'3 


UNCLASSIFIED 


,... y . .- ,; .. ... 
. .. . . 
,.:. 


UNCLASSIFIED 


self-propelled howitzers. 
Here we used values of dynamic pressure impulse 


already described. There are only eight data points. These involve 2 vehicles, 
3 orientations, 2 soil types and 7 nuclear events. Three of the data points 
(2 Badger-S and 1 Simon-7) result from shots for which there were no dynamic 
pressure measurements. Therefore, the values we used, taken from Figure 3 and 
from the scaled height of burst chart, Figure 11, are actually based upon the 
N 


displacement-inferred dynamic pressure impulse values for tanks. 


We did not attempt to use the howitzer displacement data to infer values 


of dynamic pressure impulse. To do so, we should proceed in the same manner 


as for the tank data. Referring to Figure 2, we see that of the eight data 
points only five (points 31-35) could be used. No curve of any reasonable 
degree of reliability can be drawn based on these points. 
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SECTION 3 


(This Section is Unclassified) 


RESULTS AND CONCLUSIONS 


The results described are given in Tables 2 and 3. Least squares fits 


are given in Tables 4 and 5. A description of the methods used in obtaining 
these fits and the definitions and explanations of the several measures of 
goodness of fit listed are given in the Appendix. The simplest measure of 
goodness of fit to understand is the root mLan square percent error (ED and 
Elq ). From Table 4 we see that for the displacement - dynamic pressure impulse 


fits, the errors are quite large.* 
The error in estimating dynamic pres- 


sure impulse from given displacement is fairly large (35-40%) but tolerable 
and as we have seen the values obtained, on the average, are reasonable and of 
considerable value especially when gage data are poor or lacking. 
If we were, 


however, to attempt to infer displacement from a given dynamic pressure impulse 
the errors would be much larger (150%) 
-- though again the results would repre- 


sent reasonable approximations to the average displacement to be expected. 
(The 


reason for the wide dispar A ty between the errors for the two types of inference 
is seen to be simply due to the slope of the fitted lines on the log-log plots, 
Figures I and 2, or equivalently from the exponents in the fits listed on the 
Figures and in Table 4.) 


The percent errors, E1 , shown, in Table 5 for the dynamic pressure impulse 


versus scaled ground range 
q fits are quite modest, varying from about 4 to 25%. 


These errors are quite comparable to those obtained in Reference 1 for the var- 
ious fits therein. 


The results described are shown along with the fits, in Figures I - 
10. 


These results were then used, along with the surface burst data in Reference 1 
(Figure 3.13 and in one instance Figure 3.15) to construct the desired scaled 
dynamic pressure impulse contours as shown in Figures 11 and 12. 
Also shown 


is the locus of points separating the regular and Mach reflection regions. 


We have drawn the contours to conform as accurately as possible to the 


plotted points while maintaining smoothness and a continuous variation in 
contour shape as we proceed from low to high values of scaled dynamic pressure 
impulse. 
a 


Figure 11 is the better and more reliable of the two charts. 
The curves 


fit the data points much better than for Figure 12. 
In the latter case we have 


not attempted to force the curves through the data points as this would lead to 
% 


structure in the contour shapes which is not justified by the data. 


There are various degrees of dust**, the Met-12 data corresponding to the 


heaviest dust case while the Apple 11-13 and Bee-6 data represent more moderate 


* The errors here are much larger than in Reference I because there are much 
N 


more data for trucks '(of concern in Reference 1) than for the tanks and 
howitzers of concern here. 
,*%. 


** We now believe that tle degree of "dustiness" is an indicator of the severity 


of the precursor which is related in turn to the temperature of the pre- 
shock thermal layer. 
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dust. This is the principal cause of the difficulty in Figure 12. 
However, we 


do not believe that separate charts for moderate dust and for heavy dust are 
warranted. This distinction is simply too fine considering the quality of the 
data. 


In conforming the contours to the data points in as reasonable a manner 


as possible, we note that the two charts differ for the higher scaled burst 
heights, i e., the comparable contours do not coincide, even though the same 
data are used for the two highest scaled burst heights on Figures 11 and 12. 
We cannot regard the contours in this region as well determined; our preference 
in this region is for Figure 11, because of the greater smoothness of the con- 
tours and because we believe the data at the third highest scaled burst height 
are better for Figure 11 than for Figure 12. 


Note: 
In applications of Figures 11 and 12 the scaling factor S (see 


General Nomenclature) should be used rather than W-1/3 for scaling grgund range 
and height of burst; Si should be used in conjunction with contours rather than 
W-1 /3 for instances in which ambient pressure and temperature are specified and 
differ significantly from the standard values of 14.7 psi (101.4 kPa) and 15"C, 
respectively. 


e 


16 


UNCLASSIFIED 


-I. 


I.' 


. 
- .J. 
I 


UNCLASSIFIED 


SECTION 4 


(This Section is Unclassified) 


RECOWIENDATION 


It is recomended that the Height of Burst Chart, Figure 11, be accepted 


as standard, i.e., as the most accurate obtainable from the totality of 
existent blast wave measurments. Figure 12 which includes Iq for moderate/ 
heavy dust environments is for unique conditions. The description of these 
conditions was discussed in Reference 1. Thus, Figure 12 should be used only 
for si ilar conditions or be ignored. 
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APPENDIX A 


(This Appendix is Unclassified) 


THE LEAST SQUARES PROCEDURES 


1. 
INTRODUCTION. 


All of the least squares fits in this report were obtained with the TI-59 


calculator using library programs supplied by Texas Instruments. 


All of the fits were either of the form 


y - mx + b(1 
y-ax~b 
ax 
or 
Y a a0 + a x1 + a x 2 


or ya. 1X 1 
2 2 


where y was the natural logarithm of the desired quantity; x was either the 


independent (controlled) variable or its logarithm and the least squares 


variables x1, 
x2 involved only the special case x I 
x and x2 M x2 . 


The least squares equations are obtained by minimizing the sum of the 


squares of the deviations of the data from the curve. If the experimental 


points have a variable scatter for a given small range of x, weights should 


be applied in the procedures so that the least squares equations correspond- 


ing to Equation (1), 
for example, are obtained by setting the derivatives of 


S a Mtiy 
I - 
(Mxi + b)] 
2 
(3) 


with respect to a and b equal to zero, thus obtaining two linear algebraic 


equations which can be solved for a and b. 
Here S is the weighted sum of the 


squares of the residuals about the fitted curve and 


l 
• a a2 /a i 2 


where a is a constant (to be determined from the deviations of the data from 


the fitted curve) and ai 2 is a measure of the expected deviation from the 
k 


true value for an observation y1 (of unit weight). 


In the Texas Instruments programs Wi 
1 I. In individual cases to be 


discussed we will point out that the progr=s used are nonetheless quite 


adequate for our needs. This is mainly due to tho fact that, in all cases, 


we minimized the sum of the squares of the logarithm of the desired function. 


This is equivalent to minimizing the sum or the squares of the percentage 


deviations from the curve since (tt ) - A 
. 
(4) 
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That is, for small deviations from the eux., the fractional error 
is 
y 


equal to A(bt y). 
For large deviations this is n.ot exact and in our tabula- 


tion of results we show both 


E(tn Yi - bt y. 
r 
(- 
) 


where y. 4 the value of y(xi) corresponding to the curve fit.* 


The statistical analog of Equation (4) (Reference 3) is 


Var (.n y) * 
Var y 
(Var a variance) 
. 
(4a) 


Now in all of the cases the values of y covered a large range. 
If we con- 


sider, for example, a displacement measurement of 100 metres to be in error 


by 2%, the square of its deviation from the true value is 4; if we consider a 


displacement Measurement of I metre to be in error by 20%, the square of its 


deviation from the true value is 0.04. So even though the percent errors in 


our data are not necessarily uniform (as best these errors are known), assum- 


ing them to be uniform is much closer to reality than any other assumption we 


can make. From the example just given we readily see that the coefficieots 


in the least squares fits would be determined almost entirely by the data 


with large values of y it we were to minimize absolute rather than percent 


deviations from the fitted curves. 


2. DYNAMIC PRESSURE IMPULSE VERSUS DISPLACEMENT (a) 


DISPLACEMENT VERSUS DYNAMIC PRESSURE IMPULSE (b) 


In application of our results we need fits of the data for cases (a) and 


(b), i.e., with each quantity used in the role of independent (controlled) 


and of dependent variable. 


In each subcase of this case the data are fitted very well by straight 


lines on log-log plots. Thus for fitting on the TI-S9 calculator the fit is 


of the form 
y amx b 
(1) 


with 
y -n 
q, x - tD D 
for Case (a) 


and 
y - Zn D, x - Zn lq 
for Case (b) 


and 
b=Z.nB 
. 
(s) 
- 


• In most cases we omit the subscript i denoting the ith data point; the 


summations are taken over the data points in all cases so that omission 


of the subscript i will not lead to any confusion. 
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Equation (I) can then also be written [for Case (a)], Iq - BDm 
(ia) 


with a similar relation for Case (ib). 
However, after we havo fitted Case (a) 


we can calculate (b) directly (without fitting it) provided we recall certain 


data summations from the machine memory. For Equation (1) we obtain from the 


TI-59 fit: 
rx Z -x 
Evi , 


slope a m - 
2 
(6) 


Ex. 
(rx) /N 


Yintercept " b" 
Y N' 
(7) 


(2)1/2 


and the correlation coefficient, r 
m 
(8) 


Y 
r 


r 2 is called the coefficient of determination. The various summations can be 
recalled from the machine memory. 
ax 2 and a y2 are the variances of the 


x-irray and y-array data and are given by * 
x2 
1 ~I 


a2 
12 
) 
(9) 
x 
a 
y . r(y, 
1 
(10) 


with 
oar 
t 
. 
(11) 


and 
r' 
y 
(12) 


N is the number of data points. 


Equations (6) to (8) are not symmetrical in x and y and the algebraic 


inverse of the least squares fit 
Iq- 
BDm 
(la) 


I 1/m n/ 


is 
D 
T 
but this is not a least 


squares fit 
to the data with the roles of Iq and D reversed. 


However, manipulation of Equations (6) to (10) leads to the equation 


r1 2x 
) 
1 
3 
(1x'. 
2 ]1,2 
12 


which is symmetrical in x and y, i.e., unchane 
when x and y are interchanged. 


Bar over a quantity indicates an average value. 
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Thus, if we have obtained the least squares fit Iq - BD M we know that 


the correlation ratio r has the same value with the roles of Iq and D 


reversed. It can be. shown (by taking the origin at the data centroid), 


(References 3 and 4), that the least squares fits 


Iq = BD m 
Ca) 


D = Alq " 
(ib) 


have values of Z and m which satisfy 
C m. 
(13) 


Thus, having determined m (and r) in the least squares fit to case la 


we obtain e in the least squares fit to Case (ib) from Equation (13). 
The 


constant A in Equation lb is then determined by the fact that both least 


squares fits (la) and (Ib) pass through the centroid of the data distribu- 


tion. 
(Compare Reference S.) 
Thus, from Equation (Ib) we havr 


t& D -ZJ A + f-AI q 
(lb) 


and in particular 
"-aD - Zt A C 
+ 
7 
- 
(14) 


where from Equations (11) and (12) 


Z (.eA D) 
(15) 


and these values are available from the machine memory after running case 


(1a). 
Thus from Equation (14) we determine In A, hence A and the least 


squarei fit to case (lb) is fully determined. 
-"-- 


Having obtained the least squares fits, the coefficient of determination 


provides a measure of the goodness of fit. 
By manipulation of Equation (3), 


with W. a 1, and Equations (6) tlrough (12) it can be shown that 


l1r 
2 
S 
(16) 


Noy 


In this report we are concerned only with positive correlation, so that the 


nearer r is to +1, the smaller is S, the sum of the squares of the deviations 


from the fitted curve, and hence the better the fit to the data. The expres- 


sion for 1 - 
r2 can also be put in the form 
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l- 
r 
r 
- 
(7 


thus expressing I - 
2 as a ratio of the sum of the squares of the deviations 


from the curve to the sum of the squares of the deviations from the mean (of 


the data). 
The denominator may be regarded as a normalizing factor. It is 


independent of the functional form used for the fitting function. Without 


such a denominator 1 - r2 would tend to increase with the number of data points 


- 
even if the data were excellent. 


Now letting S ., Iq and S n D be the sum of the squares of the deviations 


from the fitted curve when Zt Iq and Zn D, in turn, play the role of dependent 


variable, we have from Equation (16) 


(1- r2 ) 
Sbt Iq 
SZn D 
(18) 


N 
Na 


Here 
Z Iq =q)' 
,(...1Iq)) 


(19) 


2 
[ (Z 
D) 2 
_ 
"D) 
2 


7ZnD 
N 


in accord with Equations (9) through (12). 
Since we have seen that r is 


unchanged when the roles of D and Iq ore reversed, Equation (18) shows that 


the sums of the squares of the deviations from the curves are not the same 


for the two corresponding least squares fits (to the same data) but that 


Zn 
I 
Zn I 


Zn D 
a2 
D InID 
In all cases treated in this report we found Sbt Iq to be considerably 


less than S. a 
D" 


In an application to be made of the results of this case we wish to 


know the standard deviation of the value of the slope of the curve as well as 


the standard deviation of our observation. We use the following notation 


with the subscript c referring to the curve fit in each instance: 
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2A Iqs C 2 variance of an observation ofZn lq relative to the 


curve; (note that we have assumed the variance of 


the percent error to be independent of the abscissa 


value) 


2 
m 
a variance of slopes of curve fits. 


Then it can be shown that (Reference 6) 


a2 
S 
/N-2* 


& I qC = 
& Iq/ 


(21) 


2 
s 
/N-2 


't DC 
b D, 


S-N n 1 
N-2 
1 - 
2 a 2 


S 
- 
-. 
N 
N- 


variablez 
Iv 
the 
Iihterrlsrvrew 
loseta 


q 
q~ 


2 N-2 
a(22) 


SA I 
N-2 
I-r24nI 


Na2 
N-2 
a2 
4 
'in D 
'inD 


Using Equation(S) twice, first with D, I qbeing independent and dependent 
variable an~d then with their roles reversed, we also see that 


am 
m(3) 


m~ 


In the tabular results for the Displacement - 
Dynamic Pressure Impulse 


data curve fits, we list for each subcase the quantities N, S 
a 
I 
z q 
tAIq 


S n D' aen D,C' t, m, at. am in addition to r and the fitting functions. 


To gain further insight into the reliability of the data fits we have 


also tabulated the following additional quantities defined as follows: 


-A' 


Here there are N-2 degrees of freedom; the two degrees of freedom lost 
correspond to the number of regression coefficients (two). 
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s6 E(- 


DC 


-( 
(24) 


q 
'qC 


1/2 


D 


So 
1/2 
(25) 


I ,C a- 


and finally the root mean square percent error which is 


ED 
(S.C/N)I/2 x 100% 


(26) 


E 
- 
(S I/N) 1 / 2 x 100% 


Iq 
Iq. 


The a' quantities, Equation (25), are standard deviations of the fractional 


error (percent error apart from a factor.of 100). Thus when the number of 


data points N is fairly l-,ge (large enough so that the fit is fairly reli- 
able) the values of 
qC are almost the same as the values of ED and 


Eiq, respectively, (aside from the factor of 100 expressing the latter quan- 


tities in percent). 


The reason we have chosen to list these additional quantities, in this 


and in other data fits to be discussed shortly, is that in many cases the 


data scatter is rather large (irrespective of the functional form selected 


for fitting the data). This means that a deviation 2n y in Zt y for a given 


data point may differ considerably from Ay/y. To illustrate, suppose the 
INW 


curve fit and data point values of a displacement are 10 metres and 6 metres. 


The contribution of this point to the sum of the squares of the deviations 
of Ln D is then (tn 10- 
t 6) 2 . 0.261 (irrespective of which value is the 


curve fit value and which is the data point value). 


The contribution to S6 is 
2 


(10-6).,160P 


if the curve fit value is 10I and the data point value is 6; if, however, the 
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curve fit valut is 6 and the data point value is 10, the contribution to S 


0 6 
0.444. 


Since the least squares curve tends to pass through the data region with 


a fairly uniform distribution of points on either side, on the average the 


curve fit will be below the data point about as often as above it. 
Averaging 


th e bov 
re ult 
yi ld s 0.160 * 0.444 


the above results yields 0 
a 0.302 which is to be compared with 


the value 0.261 for the logarithmic deviation. Had we chosen the values 10 


and 9.8 instead of 10 and 6, the two compared numbers would be virtually 


identical. 


It is inherent in the nature of the data we are dealing with in this 


report that there is large data scatter while the number of data points is 


quite small. We found (Reference 1), using the same least square procedures 


described herein, that, when the scatter of the data is small, S' 1 
fairly 


close in value to S tn D and similarly Si 
is fair;r close to S 
Usually 


S < S6 nd 
SD 
< S . 
The situatlon is sometimes reversed,%owever, 
bt0 
q 
q 


usually because one or two data points with large deviations from the curve 


fit are below the curve. Conversely, for a case in which one or two data 


points are far above the curvt it can happen that S >" S t 
since there is 
-- 


no limit to the contribution such a point can make to S', y here indicating 


any quantity whose Lt has been fitted. 
(For a data point below the curve fit 


the maximum contribution to S; is 
2 
2 


Y'Yc" 
C- c) 
- 
C-.) 
= 1 .) 


The results shown in Table 4 are consistent with these statements taken from 


Refererce 1. [The large effect which a single data point can have may be 


seen by calculating 
D-D 2 
" 


- 
for point 31 of Figure 2. The value, 18.52, is 


Uc 


almost as large as S' (a 19.06) for the 8 data points of Figure 2.) 


Finally, quantities such as a' 
and Eq are directly associated with 
1q,C 


the plotted data and are easily visualized. Since they pertain to fractional 


or percent errors they apply equally well at all points along the curve fit 


(although the curve fit has greater predictive accuracy in the vicinity of 


the data centroid than toward the ends of the data range). 
On the other 
A 


26 


UNCLASSIFIED 


,...- 
, ', 
w'', 
..' ' V ' 
, 
' ",,',;, ,.',; ;., ,'-., 
.' " .' . . " .'dr. '\ ,". % -. .. .' "i 
.''. " '"'. "" -''-.-''- '.."" 
U- 


UNCLASSIFIED 


hand aft q refers to the natural logarithm of the plotted data (ordinate) 


and the latter may be quite small or large (or even negative), so that while 


ae., 
t q is a standard measure of the data scatter with respect to the fitted 


curve, it is not as simple to interpret as are the additional quantities 


which we have tabulated. 


In sum, the tabulated measures of goodness of fit provide some insight 


into the assessment of the reliability and utility of the data and the 


curve fits and help with the interpretation of results. 


The quantities described here are also tabulated for various curve fits 


of data other than displacement-dynamic pressure impulse data. We discuss 


the least squares treatment of the remaining data in the following section. 


3. ADDITIONAL LEAST SQUARES FITS. 


Scaled Dynamic Pressure Impulse versus Scaled Ground Range: 


The fitted curve is given by 


y a a0 a x + a x 2 
(27) 


where x - scaled ground range, y - b S I. 
%N. 
Sq* 


In running this case on the TI-S9 calculator we feed in the natural 


logarithm of each data point. Again, the machine generates a fit and all of 
'-' 


the relevant summtions over the data are available from the machine memory. 


For this type of fit (trivariate), however, the machine provides a quantity 
R2 rather than r as a measure of the goodness of fit. R is the multiple 


linear correlation coefficient between y and the other least squares variables 
x1 a x and x2 = x ; R2 is called the coefficient of determination (see Ref- 


erence 6) and is given by 


r 2r 2 
2r r 
I 


R2 ayl 
r 2 2 r12 
l r 2 
(28) 


- r 
. 


The subscripts 1 and 2 refer to x and x2; ry1 r, 
and r12 are given by 


equations analogous to Equation (8a). 


(,)],2-1/2 


... 


y- Y. 
" ( 
J( - 
- 
( 
(29) 
% 
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r , x 2 y - 
3xN 
] ( C r x 4 - 
2 ]XN 
( y 2 
. c N 2--] 1 / 2 


2-22 
.1/2 


r~__ 
a 
Ex3 
r.x)1)][X 
" 
:_ 
1 
. 
- 
N 


(The superficial lack of symetry between r 
and ry, results from the fact 


that the least squares variables x I 
and x2 as used here are x and x2 , respec--.,. 


tively.) The quantities ryI, ry2, and r 12 are the coefficients of correlation 
(also called simple correlation or zero-order coefficients) between y and x1 , 
y and 2. and x1 and x2. respectively. The important point here is that 


2 


( - R ) 
y 
(17a) 


R2 


so that R provides a measure of the goodness of fit precisely similar to that 


provided by r2 in the previous discussion: R2 is a measure of the closeness 


of fit of the regression plane (in It y, x 1 , x2 space)* to the data points. 


There are two further changes: replacement of N-2 by N-3 in Equations 


(21) and (2S) and replacement of I 
by S q 


it 
/N-3 


SI 
S I qC 
' 
n SjIq/ 


(21a) 


tA o,c 
n Dn 


and 


/N-3)1/2 
(2Sa) 


I'~ q c a (S~ I 


i~q 


Referring to the results listed in Table 5 we see that usually 


S. 
iI 
itq < SiIq (and in fact this relation holds for all of the cases where 


the number of data points > 11). 
The two values are usually quite close, 


however, as might be txpected, since the data dispersion is fairly modest. 


x , x • scaled ground range 
x2 a x. 
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